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ABSTRACT 

Using long-slit spectroscopy from IMoran et al.l (|2012t ) we constrain the radial depen- 
dence of the recent star formation histories of nearby galaxies with stellar masses 
greater than 10 10 M Q . By fitting stellar population models to the combination of SSFR, 
D4000 and H^a, we show that the star formation histories of many disk galaxies can- 
not be accurately represented if their star formation rates declined exponentially with 
time. Many galaxies have Balmer absorption line equivalent widths that require recent 
short-lived episodes or bursts of star formation. 

The fraction of galaxies that have experienced episodic rather than continuous 
star formation is highest for "late-type" galaxies with low stellar masses, low surface 
densities, and low concentrations. In these systems, bursts occur both in the inner and 
in the outer regions of the galaxy. The fraction of stars formed in a single burst episode 
is typically around 15% of the total stellar mass in the inner regions of the galaxy and 
around 5% of the mass in the outer regions. When we average over the population, 
we find that such bursts contribute around a half of the total mass in stars formed 
in the last 2 Gyr. In massive galaxies, bursts occur predominantly in the outer disk. 
Around a third of all massive, bulge-dominated galaxies have experienced recent star 
formation episodes that are fully confined to their outer (R > 0.7i?go) regions. The 
fraction of stars formed in a single episode is only ~ 2 — 3 % of the underlying stellar 
mass, but when we average over the population, we find that such bursts contribute 
nearly all the stellar mass formed in the last 2 Gyr. 

Recent star formation in outer disks is strongly correlated with the global atomic 
gas fraction of the galaxy, but not its global molecular gas fraction. We suggest that 
outer episodic star formation is triggered by gas accretion "events" . Episodic star 
formation in the inner regions is suppressed in galaxies with large bulge-to-disk ratio. 
This supports the idea that inner bursts are linked to instability-driven gas inflows. 

Key words: galaxy formation - optical: galaxy. 



1 INTRODUCTION 

It is now well-established that galaxies in the nearby Uni- 
verse separate rather cleanly into two classes: those with 
disky morphologies, plentiful gas and ongoing star formation 
and those that are bulge-dominated, with little gas and star 
formation, and where star formation has large l y ceased (e.g. 
IStrateva et~aH 1200 ll : iKauffmann et all l2003bl : iBaldrv et all 
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|2004 ). Most of the baryons in the Universe are not locked up 
in stars, but reside in the circum-galactic medium associated 
with virialized dark matter hal os or in a more diffuse phase 
in th e intergalactic medium (|Fukugita. Hogan. fc Peebles! 
Il998h . It is thus reasonable to postulate that the observed 
bimodality in galaxy properties is in some way linked with 
the ability of galaxies to accrete gas from the external envi- 
ronment. 

It is very difficult to probe gas infall onto galaxies di- 
rectly. Resolved studies of the stellar populations of galax- 
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ies may, however, provide som e important clu e s as t o how 
accretion may refuel galaxies. iBell fc de Jond (|2000l ) used 
a large sample of low inclination spiral galaxies with radi- 
ally resolved B,R,K photometry to derive estimates of age 
and metallicit y as a function of position within the galaxy. 
IBell fc Bowel (|200d ) then interpreted these results using 
very simple models of disk galaxy evolution. They showed 
that a "closed box" model in which the star formation his- 
tory and the metallicity of a given area in any galaxy de- 
pend only on the initial local gas surface density did not 
provide a good fit to the observational results. Real galax- 
ies had stronger age gradients than predicted, and the pre- 
dicted slope of the age/central surface brightness correlation 
was steeper than observed. Models which included infall of 
gas (with an infall rate that was larger for low mass galax- 
ies) were required to fit the age/surface brightness relations. 
Outflows were required to explain the observed correlation 
between stellar metallicity and galaxy mass. 

IMunoz-Mateos et al.l (|2007T ) examined the specific star 
formation rate profiles of a sample of nearby, face-on spi- 
ral galaxies with UV ph otometry from the GA LEX At- 
las of Nearby Galaxies (|Gil de Paz et al.l 120071 ) and K- 
band photometry fro m the Two Micron All Sky Survey 
(|Skrutskie et alj|2006h . They showed that on average, galax- 
ies were forming stars at higher relative rates in their outer 
regions compared to their inner regions. More recently, 
IWang et al.l (|201ll ) showed that UV/optical colour gradi- 
ents of galaxies are strongly correlated with their atomic 
gas content. Gas-rich galaxies with high Hi gas masses had 
stronger colour gradients than galaxies of the same stellar 
mass, size, and NUV- r colou r wit h average Hi c onten t. Both 
IMunoz-Mateos et al.l |2007h and I Wang et all (|201ll ) inter- 
preted their results in the context of an "inside-out" picture 
of disk galaxy formation, which has commonly served as a 
basis for semi-analytic models of the formation of discs in 
the context of cold dark matter cosmologies. 

The UV/optical colours of galaxies depend strongly on 
the amount of dust extinction in the galaxy and also are 
sensi tive to star forma tion occurring over timescales of ~1 
Gyr (iKong et al.l l2004f). They thus cannot be used to assess 
whether the recent star formation history of a galaxy has 
been smooth or episodic. Combinations of narrow-band stel- 
lar absorption line indices, in particular the 4000 A break 
strength and the equivalent widths of Balmer absorption 
lines , such as H7 or H<5, are able to probe the recent 
star formation histories of galaxies in more detail. Galax- 
ies that have undergone a burst of star formation in the last 
1-2 Gyr will have stronger Balmer absorption line equiv- 
alent widths for a given value of 4000 A break strength, 
compared to galaxies that have been forming their stars 
continuously jDressler fc Gunnlll983l ; IPoggianti et al.|[l999l ; 
iKauffmann et alJl2003al ). 

This paper aims to constrain the radial dependence of 
the recent star formation histories of disk galaxies. We com- 
pare spectral indices, the 4000A break strength (D n 4000), 
a Balmer line index derived using an approach based on 
principal component analysis (PCA), and the present-day 
star formation rate over stellar mass (specific star forma- 
tion nite_j_SSFll)_toa librar y of models generated from 
the iBruzual fc Charlotl (|2003l ) population synthesis code. 
We use the best fit models to constrain the timescale over 
which stars have been formed at different radii in the disk. 



We show that the star formation histories of many spiral 
galaxies cannot be be described by simple exponentially de- 
clining star formation histories. Additional recent (< 2 Gyr 
old) episodes of star formation are required in order to fit 
D n 4000, H7A and SSFR simultaneously. The mean fraction 
of stars formed in these recent episodes is similar in the in- 
ner and outer regions of galaxies, except for massive galaxies, 
where the bursts occur predominantly in the outskirts of the 
galaxy. 

We also study how the recent star formation histories 
and gas phase metallicities in the inner and outer regions of 
the galaxies in our sample correlate with their atomic and 
molecular gas content. Our main result is that the recent star 
formation history and the metallicity in both the inner and 
in the outer regions of galaxies are strongly correlated with 
their total atomic gas mass fraction. In contrast, only the 
inner star formation history and metallicity are correlated 
with total molecular gas fraction. 

Our paper is organized as follows. Section 2 intro- 
duces the 21cm and C O(l-O) line data we use from GASS 
dCatinella et alj|20ich and COLD GASS |Saintonge et al.l 
120111 ) surveys and the long-slit spectr oscopic observation s 
from the MMT follow-up observations (|Moran et al.ll2012T l. 
We explain how we create a library of model spectra and use 
it to extract estimates of parameters such as burst mass frac- 
tion in Section 3. Based on the results of the fits, we classify 
our sample galaxies into five groups according whether their 
star formation histories are best described by continuous or 
burst models in their inner and outer regions. In Section 4, 
we examine how the fraction of galaxies in these different 
classes depends on their location in the plane of stellar mass 
versus stellar surface density/concentration/colour. We also 
examine how the fraction of stars formed both continuously 
and in bursts in the inner and outer regions of the galaxy 
depends on the location of galaxies in these planes. We sum- 
marize and discuss our results in Section 5. 



2 DATA 

The data set on whi ch this analysis is based has been de- 
scribed in detail in (jMoran et alJl2012h . The original par- 
ent sample consists of galaxies with Hi line flux data from 
the Arecibo telescope that were observed as part of the 
GALEX Areceibo SDSS (GASS) survey (Catinella et al. 
2010; 2012). The reader is referred to these papers for de- 
tails about the sample selection and the Hi observations. 
A subset of 300 galaxies were selected from this sample for 
CO ( J = 1 — 0) observ ations using the IRAM 30m telescope 
(|Saintonge et al.ll201ll ). Of these, around 200 galaxies now 
have reduced long-slit spectroscopy obtained from both the 
Blue Channel Spectrograph on the 6.5m MMT telescope on 
Mt. Hopkins, AZ, and the Dual Imaging Spectrograph on 
the 3.5m telescope at Apache Point Observatory (APO; see 
iMoran et al. I l2010l for details about the observational setup 
and data reduction). The wavelength ranges for MMT and 
APO observations are ~3900-7000 A at a spectral resolu- 
tion of ~4 A (90 km s _1 ), and ~3800-9000 A at a spectral 
resolution of ~6-8 A (150 km s _1 ). 

The spectra were spatially binned outward from the 
galaxy center to ensure an adequate S/N in each spatial 
bin. All bins have minimum extent equal to the slit width: 
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1.25" for the MMT observations) and 1.5" for the APO ob- 
servations. The bin size is typically <3" , corresponding to 
a physical size of ~1.5-3 kpc at the redshift of the galax- 
ies in our sample. The S/N per bin is required to be larger 
than 15 (per A) in the inner regions of the galaxy. In the 
outer, lower surface brightness regions, these requirements 
are progressively loosened in order to compromise between 
S/N and spatial resolution: we adopt S/N >10 for bins out 
to a distance of 4.5", S/N >8 for bins between 4.5 and 6" 
and finally S/N >6 in the outermost bins. The spectra were 
flux-calibrated by matching to SDSS g and r-band photom- 
etry measured through an aperture matched to the slit. 

All galaxies have NUV and FUV fluxes from the 
GALEX medium imaging survey. Global galaxy parameters 
such as stellar mass, stellar surface mass density, and concen- 
tration index were taken from the MPA/JHU value-added 
catalogs ( |http: / /www.mpa-garching.mpg.de/SDSS| ) . 



3 ANALYSES 

3.1 A PCA-based Balmer absorption line index 

The use of "classic" Balmer absorption line Lick indices such 
as H5.4 or H7A measured directly from the galaxy spectra 
can be problematic for two reasons: 1) Because these indices 
are defined over a narrow range in wavelength, the signal-to- 
noise of the measurements is quite low, especially towards 
the outer regions of the galaxy, where the spectra may only 
have S/N ~ 2 per pixel. 2) In star-forming galaxies, the 
Balmer absorption lines are filled in by emission. A robust 
measurement of the absorption line equivalent width hinges 
on subtracting this emission accurately. 

If stellar population models fit the observed spectra ac- 
curately, they provide one way of solving the problem. Be- 
cause the models are fit to the spectra over a wide range in 
wavelength, the errors on Balmer line indices derived from 
these fits dec rease with respec t to those for directly mea- 
sured indices (jWild et al.ll2007l ) . In addition, the wavelength 
regions where emission may be present are masked when the 
fitting is carried out. 

In this paper, we apply the principal compon ent analy- 
sis (PCA) method described in lChen et al.l ll2012h. which i s 
an extension of the method described in IWild et al.l (|2007h . 
The principle of PCA method is that a galaxy spectrum 
can be decomposed into a set of orthogonal principal com- 
ponents (PCs), which can be linked to physical properties 
such as fraction of stars formed over the last Gyr, or to es- 
timates of indices such as D4000, H<5a, H7A using a library 
of star formation histories generated using population syn- 
thesis models. By comparing the amplitudes of the principal 
components of the observed spectra to those of the modelled 
spectra, the likelihood distribution of such parameters can 
be calculated and thus constrained. The advantage of this 
approach is that it improves our estimates of Balmer absorp- 
tion line equivalent widths, particularly for the spectral bins 
with poor S/ N in the ou t skirts of the galaxies. The reader 
is referred to IChen et al.1 (|2012T ) for further details. In this 
work, we find the linear combination of PC amplitudes that 
best represents the indices D4000 and H7A for each spec- 
trum. A probability distribution function (PDF) for either 
D4000 or H7A can be built using \ 2 as weights. We use the 



median value of the PDF as our nominal estimate of the pa- 
rameters, and the 16th to 84th percentile range of the PDF 
as the ±1(7 confidence interval. 

3.2 Comparison with a Library of SEDs of Model 
Galaxies 

We compare two PCA-derived spectral indices, which we will 
denote PCA Dn ( 4 ooo) an d PCAh 7a , as well as the present-day 
star formation rate over stellar mass (SSFR) derived from 
the ex tinction-corrected Ha flux as described in lMoran et al.l 
(|2012T l , to a library of models generated from the Bruzual & 
Chariot (2003) population synthesis code. As we will demon- 
strate, these three quantities jointly place strong constraints 
on whether the galaxy has been forming stars smoothly or 
episodically over the past 1-2 Gyr. 

We create a l ibrary of models by using the population 
synthesis code of iBruzual fc Chariot! f2003). The underly- 
ing model has a continuous SFR, declining exponentially 
with time SFR(i) oc exp(-ft), with 7 uniformly distributed 
between (i.e. constant star formation rate) and 1 Gyr -1 . 
Stars begin to form at a lookback times between 13 Gyr and 
1.5 Gyr. 

We then superimpose additional bursts of star forma- 
tions onto these continuous models. The strength of the 
burst is defined as the fraction of the stellar mass produced 
by bursts relative to the total mass formed by continuous 
models, and is distributed logarithmically between 0.001 and 
4. Bursts are described by top-hat functions and have du- 
rations between 3xl0 7 to 3x10 s years. The combination 
of Balmer absorption lines and 4000 A break strength is 
not sensitive to star formation episodes that have occurred 
longer than about 2 Gyr ago, so the bursts occur randomly 
between lookback times of 2 Gyr and the present. The metal- 
licity range of our models is distributed uniformly from 0.005 
to 2.5Zq. Each model inc ludes dust extinction b ased on the 
two-component model of ICharlot fc Fan] (|2000l ). where the 
V-band optical depth t v follows a Gaussian distribution with 
a peak at 1.78 and a width of a = 0.55. The t v distribution 
is truncated so that it only spans the range 0-4. fi is the 
fraction of that optical depth affecting stellar populations 
older than 10 Myr and is uniformly distributed between 
and 1. 

3.3 Fitting procedure 

We first discard those spectral bins with contamination from 
background or foreground sources (3 spectral bins in total, 
originating from 2 galaxies). We also exclude the spectral 
bins with large errors in the PCA indices — PCAh 7a > 1 
and PCA Dn ( 4 QQQ)> 0.2. This leaves a total of 2725 spectral 
bins out of the original 2812 bins. We first check whether the 
2 PCA indices and the SSFR can be fit using the library of 
continuous SF models. The bins where such fits fail are then 
fit with the full library of continuous and burst models. We 
have used both the best-fit model with minimum chi-square 
to estimate quantities such as the stellar mass formed in 
the past 2 Gyr, as well as the median of the probability 
distribution function. In practice, the two estimates yield 
virtually identical results (see Figs 1 and 2). We record the 
stellar mass produced in the last 2 Gyrs. When the continu- 
ous models do not provide acceptable fits within la error of 
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Figure 1. Model fits for galaxies with star formation histories of Type A (columns 1 and 2) and Type B (columns 3 and 4). The first 
row shows SDSS cut-out images of the galaxy. Rows 2,3, and 4 show the radial run of the specific star formation rate, the 4000 A break 
index and the ^[7^4 index. Black points show the data along with error bars. Red/blue lines are the model fits using minimum chi-squarc 
model and the median of the PDF, respectively. In the bottom row we plot the look-back time when half of stellar mass was formed. 
Note that radii are scaled by dividing by R90, the radius enclosing 90% of the r-band light. 
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Figure 2. Examples of a galaxy of Type C (Columns 1-2) and Type D ( Columns 3-4). The rows and line styles have the same meaning 
as in Figure 1. Columns 1 and 3 show the best-fit continuous models. Columns 2 and 4 show the best-fit models that include bursts. 
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data, we also record the stellar mass produced in the burst 
mode during this time period. 

3.4 Separating inner and outer regions of galaxies 

In the following sections, we will often refer to the mass of 
recently formed stars in both the "inner" and the "outer" 
regions of the galaxy. We choose 0.7Rgo as the nominal ra- 
dius to partition our galaxies into inner versus outer spec- 
tral bins, where Rgo is the radius enclosing 90 percent of 
the r-band Petrosian flux. We weight each spectral bin by 
the stellar mass enclosed by the circular annulus formed by 
the inner and outer boundaries of the bin. We then calcu- 
late weighted averages of quantities such as the fraction of 
recently-formed stellar mass and gas-phase metallicity. 

3.5 Categorization of Galaxies According to Their 
Star Formation Histories 

We divide our galaxy sample into five classes according to 
their star formation histories in their inner and outer re- 
gions: 

(1) Type A: Those that are well described by continuous 
models at all radii, i.e. the star formation histories are well 
described by a varying exponential decline time as a function 
of radius. Type A galaxies have t s tart ranging between 13 
and 12 Gyr in lookback time. 

(2) Type B: Those that are well described by continuous 
models at all radii, but with tstart varying between 12 and 
2 Gyr. 

(3) Type C: Those that are well described by contin- 
uous models with tstart varying between 13 and 2 Gyr in 
their inner regions, but which require a burst in their outer 
regions. 

(4) Type D: Those that require bursts in both their inner 
and their outer regions. 

(5) Type E: Those that require bursts only in their inner 
regions. 

The number of galaxies of type A, B, C, D, and E is 45, 
24, 34, 104, 2 respectively. Some examples of each type are 
presented in Figs [T] and [5] Fig. Q] shows examples of galaxies 
of Types A and B, which are well fit by pure exponential 
models. In the first row, we show SDSS cut-out images. The 
2nd, 3rd, 4th rows show SSFR, D4000 and H7A as a function 
of radius in the galaxy. Note that the radii have been scaled 
by dividing by Rgo. In these panels, black points show the 
data along with the errors. Red/blue lines are model fits 
using the minimum chi-square models and the median of the 
PDF, respectively. In the 5th row, we plot the parameters 
of the minimum chi-square models as a function of radius. 
The 5th row shows the lookback time when half of stellar 
mass was formed at each radius. 

Fig. [2] shows examples of galaxies of Types C and D, 
which require bursts. The rows are the same as in the pre- 
vious figure. In the left column, we show the best-fit con- 
tinuous models. As can be seen, in both cases such mod- 
els generally work better near the centers of galaxies. In 
the outer regions, the deviations in the fit to PCAh 7a and 
PCA Dn (4ooo) become large. The right panels show how addi- 
tion of a burst bring the models into much better agreement 
with the data. We note that the physical size of each bin for 
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Figure 3. The average deviation of the observed values of 
D n (4000), H7A and SSFR from the values predicted by the best- 
fit continuous models for all Type C (upper panel) and Type D 
(lower panel) galaxies as a function of radius. Note that the devi- 
ations are scaled by the measurement errors a and are plotted in 
units of a. This takes into account the fact that the S/N of the 
binned data is higher in the inner regions than the outer regions. 



our sample is at least ~1.5-3kpc, and the discrepancies be- 
tween the data and the continuous models often span more 
than 2 bins. This means that the bursts are not generated 
by individual Hn regions, but span a significant portion of 
the outer disk. 

In Fig. we plot the average deviation of the continu- 
ous models from the observations as a function of radius for 
Type C (upper panel) and Type D (lower panel) galaxies. 
In order to account for the fact that the observational errors 
increase as function of radius, we scale the deviation by a 
before taking the average . As expected, the centers of Type 
C galaxies are well fit, with the average deviation increasing 
towards the outer regions. The lower panel of Fig. [3] shows 
that the same general trend is also true for Type D galaxies. 
In the central regions, the average deviation from a continu- 
ous model is only ~ la, but in the outer regions the average 
deviation increases by a factor of 2-3. As we will show in 
the next section, the main difference between galaxies with 
Type C and Type D star formation histories is stellar mass. 
Galaxies with recent star formation only in their outer re- 
gions are predominantly systems with stellar masses greater 
than ~ 10 n M Q . 
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Figure 4. Left: Fraction of galaxies with different types of star formation history in the two-dimensional plane of stellar surface density 
versus stellar mass. Poisson error bars are shown. Right: Fraction of the stellar mass formed in the last 2 Gyr in the inner region 
(R < 0.7-Rgo, red), outer region ( R > 0.7i?9o, cyan), and in the whole galaxy (green). The hatched area shows the fraction of stellar 
mass formed in the last 2 Gyr in the burst mode. Error bars have been computed using boot-strap resampling. 
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Figure 5. As in FigU except for galaxies in the two dimensional plane of concentration versus stellar mass. 



4 RESULTS ratio of the galaxy. Finally, the global colour is a measure of 

the average age of the stellar population. 

4.1 Dependence of star formation history on 
stellar mass, surface density, concentration 
and colour 

In this section, we investigate how the fraction of galax- 
ies with different types of star formation histories defined In order to understand whether the star formation his- 

in Section 3.5 depends on stellar mass, galaxy colour and tory of a galaxy is determined primarily by its mass or by its 

structural parameters such as surface density and concen- structural properties, we will analyze SFH trends in two di- 

tration. Stellar mass is a good tracer of the mass of dark mensional planes of stellar mass versus surface density and 

matter halo in which the galaxy resides. The surface density concentration. We will also examine how the fraction of stars 

of the galaxy measures the degree to which the gas has lost formed both continuously and in bursts in the inner and 

angular momentum prior to being converted into stars. The outer regions of the galaxy depends on the location of the 

concentration index is a good proxy for the bulge-to-disk galaxies in these 2D planes. 
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Figure 6. As in Fig. [4] except for galaxies in the two dimensional plane of NUV-r colour versus stellar mass. 



4-1.1 The plane of stellar mass versus stellar surface 
density 

All galaxies are divided into 6 bins by dividing their stellar 
mass into 3 bins, 10.0-10.4 M , 10.4-10.8 Mq, and 10.8- 
11.2 Mq in logM*, and surface mass density into 2 bins: 
log/x* < 9 Mq kpc -2 and log/i, > 9 Mq kpc~ 2 . In each 
bin, we calculate (1) the fraction of galaxies with different 
types of star formation histories (left panel ), (2) the fraction 
of the stellar mass formed within the last 2 Gyr in the inner 
region, the outer region and in the entire galaxies (right 
panel), (3) the fraction of the stellar mass formed in the 
burst mode in the last 2 Gyr in the inner region, the outer 
region, and in the entire galaxy (hashed histograms in the 
right panel). 

As shown in the left panel of Fig. [4] the star formation 
histories of galaxies clearly depend on both stellar mass and 
surface density. The fraction of quiescent galaxies without 
bursts (types A and B) depends mainly on surface density 
rather than stellar mass. There is some tendency for quies- 
cent galaxies to be somewhat more active (i.e. be of Type 
B) if they have low masses. 

The fraction of galaxies that have experienced recent 
bursts in both their inner and their outer regions (type D) 
decreases both as a function of stellar mass and stellar sur- 
face density. In the stellar mass bin with log M, in the range 
10-10.4, type Ds make up more than two thirds of the entire 
population, but type Ds make up only a small fraction of the 
galaxies our highest M*/^* bin. Interestingly, however, the 
fraction of galaxies that have experienced recent bursts in 
their outer regions only (Type Cs) exhibits the opposite de- 
pendenc of on the mass of the galaxy. In the highest stellar 
mass bin with logM* in the range 10.8-11.2, there are as 
many Type C galaxies as there are quiescent (Type A and 
B) galaxies. 

The right panel of Fig. [4] displays the fraction of stars 
formed in the last 2 Gyr as a function of location in the 
plane of stellar mass and stellar surface density. We show 
the fraction of recently formed stars in the whole galaxy (in 
green) as well as the fraction in the inner (R < 0.7i?go, red) 



and outer (R > 0.7i?go,blue) regions. The fraction of stars 
formed in the last 2 Gyr in the burst mode is indicated as 
hashed shading. As expected, the fraction of recently formed 
stars is highest for galaxies with low masses and low surface 
densities and lowest for galaxies with high masses and high 
surface densities. Interestingly, the recent star formation is 
always rather evenly split between the inner and outer parts 
of galaxies, except for the most massive galaxies with low 
stellar surface densities, where more star formation occurs 
in the galaxy outs kirts. These r esults are consistent with a 
recent analysis by I Perez et al.l (|2013T ) of star formation in 
the inner and outer regions of galaxies using a sample of 
~ 100 galaxies with integral field unit (IFU) spectroscopy. 
What is new here, is that we can estimate the fraction of 
the recent star formation that occurred in bursts. We see 
that the burst fraction is highest in massive galaxies. The 
majority of galaxies with stellar masses in the range 10-10.4 
in log M 4 require bursts in both their inner and in their outer 
regions (i.e. they are classified as type D systems), but the 
total mass of stars formed in the bursts over the last 2 Gyr 
is comparable or smaller (in the case of low mass galaxies 
with low densities) than that formed in the quiescent mode. 
In high mass galaxies, the opposite is true - the majority of 
the recent star formation is contributed by bursts. 



4-1.2 The plane of stellar mass versus concentration 

Fig. [5] is similar to Fig. 2] except that we divide galaxies 
according to concentration index (Rgo/Rso) at a value of 
2.6. As can be seen, the results are quite similar to those 
obtained when galaxies are divided by stellar surface mass 
density. One notable difference is that in Fig. [5] we see 
that all type A have the division between quiescent galaxies 
and galaxies with actively ongoing star formation is much 
more pronounced. We see that all Type A galaxies have 
Rqq/Rzm > 2.6, whereas Fig. [4] shows that there are a 
signficant number of Type A galaxies with stellar surface 
mass densities less than 10 9 Af© kpc -2 . We note that al- 
though high concentration index (equivalently bulge-to-disk 
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ratio) rather than high stellar surface density (equivalently 
contraction factor) appears to be the necessary condition 
for galaxies with quiescent Type A star formation histo- 
ries. Although high bulge-to-disk ratio is a necessary condi- 
tion, it is not a sufficient condition - a sizeable fraction of 
bulge-dominated galaxies have experienced bursts of star- 
formation in the last 2 Gyr. 

4-1.3 The plane of stellar mass versus colour 

In Fig. |6]we divide galaxies according to (NUV-r) colour at 
a value of 4.8 . As can be seen from the right-hand panel, 
red global colours select out the galaxy population that have 
experienced very little recent star formation both in their in- 
ner and in their outer regions. It is noteworthy that almost 
all galaxies with quiescent star formation histories (Types 
A and B) belong to the red population with NUV-r > 4.8. 
Galaxies that have experienced bursts in their outer regions 
are roughly split evenly between the red and the blue popu- 
lations, while type D galaxies that have experienced bursts 
in both their inner and outer regions are predominantly blue. 
We note, however, that a significant population of red Type 
D galaxies is found in the lowest stellar mass bin. 

4-1-4 Distribution of SFH parameters of galaxies 

In the previous sub-sections we studied the fraction of galax- 
ies with bursty versus continuous star formation histories as 
a function of stellar mass, colour and structural parameters. 
We also looked at the contribution of bursts to the recent 
star formation in the inner and the outer regions of galaxies, 
as a function of these same parameters. In this sub-section, 
we analyze the distribution of burst strengths in galaxies, as 
well as the look-back time when half the stars were formed. 
For simplicity, we confine our attention to the two dimen- 
sional plane of concentration versus stellar mass. 

(i) The Fraction of stars formed in bursts. Fig. [7] shows 
the distribution of the fraction of stars formed in bursts in 
the inner and outer regions of galaxies as a function of po- 
sition within this plane . As seen in the top panel, the frac- 
tions of stars formed in bursts in the inner region is generally 
around a few percent for bulgy galaxies (-R90/-R50 > 2.6) or 
for massive galaxies (logAf* > 10.8). In the inner regions 
of less massive, disk-dominated galaxies the median burst 
strength is ~ 10% and the distribution is much wider. The 
bottom panel of Fig. [7] shows that burst strengths in the 
outer regions of galaxies are insensitive to whether it is bulge 
or disk-dominated. The main controlling factor is the mass 
of the galaxy. In the highest stellar mass bin, the median 
burst strength in the outer region is ~ 3 %, about a factor 
of two lower than in the other two bins. We also note that 
in massive galaxies, the outer bursts are stronger than the 
inner bursts on average, whereas the opposite is true for low 
mass galaxies. 

(ii) The look-back time when half the stellar mass was 
formed. In Fig. [5] we make similar plots of the look-back 
time when half the stellar mass was formed (T50), which pro- 
vides information about the mean stellar ages of the galax- 
ies. The trends in this plot are much weaker. Only a very 
small fraction of galaxies have T50 less than ~4x 10 9 Gyr 
in either their inner or their outer regions. This is in keeping 



with the results in the previous section, which shows that 
although recent bursts are required to fit the spectral param- 
eters, the fraction of mass formed in these bursts generally 
does not exceed ~ 10 — 20% of the total stellar mass. Disk- 
dominated galaxies generally have younger mean stellar ages 
than bulge-dominated galaxies of the same stellar mass, and 
the most massive galaxies have the oldest mean stellar ages. 
The trend for the most massive, bulge-dominated galaxies 
to have the oldest mean stellar ages holds in both the inner 
and outer regions of the galaxies, but is somewhat stronger 
in the outer region. 

We also note that outer star formation in massive galax- 
ies discussed in the previous section, is manifested in this 
plot as shift in the lower 25th percentile of the mean age 
distribution, rather than a shift in the median value. This 
means that outer disk growth is a phenomenon that per- 
tains to a minority of the massive galaxies, rather than to 
the population in general. 



4.2 Relations between star formation history, 

HI/H2 mass fraction and gas-phase metallicity 

In the previous section, we explored how the star formation 
histories in the inner and outer regions of galaxies depend 
on their stellar mass, structural properties and colour. In 
this section, we examine how star formation history depends 
on the gas con tent and gas-phase metallicity of a galaxy. In 
previous work. lMoran et al.l (|2012h showed that around 10% 
of all galaxies with stellar masses greater than log A/* > 10 
exhibit strong drops in gas-phase metallicity in their outer 
regions. These galaxies tended to have high Hi content and 
to have actively star-forming outer disks. 

In this analysis, we explore the correlations between the 
fraction of recently-formed stars and gas-phase metallicity 
in the inner and outer regions of the galaxy with both the 
Hi and H2 gas mass fraction. We also explore how the cor- 
relations change as we vary the radius R cu t used to divide 
the inner galaxy from the outer galaxy. 

The spectral bins with only upper or lower limits on 
metallicity are discarded before we carry out our analysis. 
For each value of R cu t (in units of R90), we compute corre- 
lation coefficients of Spearman rank tests for 

• The relations between the fraction of stars formed in 
the inner/outer regions of the galaxy and total Hi gas mass 
fraction. 

• The relations between the fraction of stars formed in 
the inner/outer regions of the galaxy and total H2 gas mass 
fraction. 

We only take into account the correlation coefficients with 
probability larger than 0.97. For the coefficients with prob- 
ability smaller than 0.97, we set the coefficients to zero. (In 
practice, this only happens for the relations between outer 
star formation and metallicity and H2 gas mass fraction 
when R C ut is greater than Rgo-) 

Fig[5] summarizes the results of this exercise, rcut is in 
units of R/R90. The upper and lower panels show the results 
for the correlations with Hi and H2 gas mass fraction, re- 
spectively. Coefficients with values larger and smaller than 
zero indicate positive and anti-correlations between the pa- 
rameters. 
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Figure 7. Top: Distribution of the fraction of stars formed in bursts in the inner regions of galaxies in the two-dimensional plane of 
concentration versus stellar mass. Bottom: As in the top panel, but in the outer region. The red, blue, green dotted lines indicate 25%, 
median, 75% values of the distribution in each panel. 



As can be seen in the upper panel of the figure, the frac- 
tion of recently-formed stars and the gas-phase metallicity 
in both the inner and the outer regions of the galaxy corre- 
late strongly with the global atomic gas mass fraction. The 
strength of the correlations does not depend on the adopted 
value of R cu t, until R cu t reaches values greater than 0.9 R90. 
At this point, stronger correlations between the gas-phase 



metallicity and the fraction of stars formed in the outer re- 
gions and the total Hi content are seen. Interestingly, the 
correlation between the fraction of stars formed in the inner 
region of the galaxy and total Hi gas fraction is independent 
of the value of R cu t- Even as R cu t decreases to values near 
zero, there is still a correlation between the fraction of stars 
formed in the very inner regions of the galaxy and its total 
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Figure 8. Top: Distribution of the look-back time when half the stellar mass was formed in the inner regions of galaxies in the two- 
dimensional plane of concentration versus stellar mass. Bottom: As in the top panel, but in the outer region. The red, blue, green dotted 
lines indicate 25%, median, 75% values of the distribution in each panel. 



atomic gas content. We will discuss possible reasons for this 
in the next section. 

In the lower panel of Fig. [9] we see that the gas-phase 
metallicity and the fraction of recently-formed stars in the 
inner regions of the galaxy are also both correlated with H2 
mass fraction. The correlation of these inner quantities with 
H2 mass fraction is stronger than the correlation with Hi 



mass fraction (r values of ~ 0.6 rather than ~ 0.4). Similar 
to what is seen for the atomic gas, the strength of the cor- 
relations of the inner quantities with H2 gas mass fraction 
are independent of R cu t- The outer quantities, however, are 
much more weakly correlated with H2 mass fraction and for 
Rent > 0.7, the correlations with H2 mass fraction disappear 
entirely. 
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Figure 9. Top: The correlation coefficients of Spearman rank 
tests for the fraction of recently-formed stellar mass and gas-phase 
metallicity as a function of HI mass fraction. The four parame- 
ters, f2gjn, f2g_ot, Z_in, Z_ot, are M^Gj/r I 'M ,„„ for the inner 
region, yi-„2Gyr/^*now for the outer region, inner metallicity, 
and outer metallicity. The correlation coefficients are plotted as 
a function of R cu t, the radius (in units of Rgo) used to separate 
the inner region of the galaxy from the outer region. Bottom: As 
in the top panel, but for the four parameters as a function of H2 
mass fraction. 



5 SUMMARY AND DISCUSSION 

We now summarize our main results as follows: 

(i) By fitting stellar population models to the combina- 
tion of SSFR, D4000, H7A, we show that the star formation 
histories of many disk galaxies cannot be accurately repre- 
sented if the star formation rate has declined exponentially 
as a function of time. Many galaxies have strong Balmer ab- 
sorption lines that require recent short-lived episodes of star 
formation in both their inner and in their outer regions. 

(ii) The fraction of galaxies that have experienced such 
episodes in both their inner and outer regions is highest for 
systems with low stellar masses, low surface densities, and 
low concentrations, i.e., late-type galaxies. 

(hi) Around a third of all massive (logAf, > 10.8), 
bulge-dominated galaxies have experienced recent star for- 
mation episides only in their outer regions. 

(iv) For low mass, disk-dominated galaxies, the fraction 



of stars formed in a single burst episode is typically around 
15% of the total stellar mass in the inner regions of the 
galaxy and around 5% in the outer regions of the galaxy. 
When we average over the population, we find that such 
bursts contribute around a half of total mass in stars formed 
in the last 2 Gyr. 

(v) For massive galaxies (logM, > 10.8), the fraction 
of stars formed in a single episode is only ~ 2 — 3 % . Aver- 
aging over the population, however, we find that such bursts 
contributed nearly all the mass in stars formed in the last 2 
Gyr. 

(vi) The amount of recent star formation in both the 
inner and outer regions of a galaxy is positively-correlated 
with its total atomic gas content. In contrast, only the inner 
star formation is correlated with total molecular gas content. 

(vii) Very similar results are obtained for gas-phase 
metallicity. Metallicity in both the inner and outer regions 
of galaxies are negatively correlated with global atomic gas 
fractions. The metallicity in the inner region of the galaxy 
is negatively correlated with total the molecular gas mass 
fraction. 

We hypothesize that these results can be understood 
if galaxies accrete atomic gas and form stars episodically. 
The fact that episodic component of the star formation oc- 
curs primarily in the outer regions of massive galaxies (right 
panels of Figs. 0] and O argues that gas may be accreted 
at large radii. We note that the non-episodic (continuous) 
component of the star formation is much more evenly split 
between the inner and the outer regions in these systems. 

In low mass galaxies, both the episodic and the continu- 
ous components of the star formation are evenly divided be- 
tween the inner and outer regions of the galaxy. It is tempt- 
ing to postulate that this may indicate that gas accretion 
occurs in a different mode in low mass systems - some the- 
oretical models predict that gas should be accreting in the 
form of "cold streams" in dark matter halos with low mass, 
and from quasi-static a tmospheres of hot gas in dark matter 
halos of high mass (e.g. lKeres et al.ll2005l ; lDekel fc Birnboiml 
l2006h . 

However, our data does not allow us to rule out the 
possibility that the inner and outer bursts in galaxies may 
have different triggers, for example gas accretion may be the 
trigger in the outer regions and disk-driven instabilities in 
the inner regions. Indeed, at fixed stellar mass, the fact that 
burst strengths are larger in the interiors of disk-dominated 
galaxies densities (Fig. top panel), argues that disk in- 
stabilities may indeed be key to the origin of the episodic 
nature of the inner star formation. The Toomre Q parameter 
scales inversely with disk surface density. 

The idea that gas accretion at large radii is followed by 
instability-triggered radial inflows of gas and a burst of star 
formation in the central regions of the galaxy is also consis- 
tent with our result that gas-phase metallicity and the frac- 
tion of recently formed stars in the inner region of the galaxy 
is best correlated with its molecular gas content. Molecular 
gas extends to much smaller radii than the atomic gas in all 
disk gala xies and also exist s at considerably higher surface 
densities (|Lerov et alj|200"3 ). Why then do metallicity and 
recent star formation in the inner regions of galaxies also cor- 
relate (albeit more weakly) with the total atomic gas content 
of the galaxy? The most reasonable explanation, in our view, 
is that the atomic gas represents a longer-term reservoir for 
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star formation in the inner galaxy. In other words, if the 
galaxy has no atomic gas, there is simply no raw material to 
be transported inwards to be converted into moelcular gas 
and stars at galactic centers. 

Finally, we would like to note that we began this paper 
with the statement "It is now well established that galax- 
ies in the nearby Universe separate rather cleanly into two 
classes: those with disky morphologies, plentiful gas and on- 
going star formation and those that are bulge-dominated, 
with little gas and star formation, and where star formation 
has largely ceased." Examination of the left panel of Fig. [4] 
indicates that even at high stellar masses, bulge-dominated 
galaxies have actually experienced a wide variety of differ- 
ent star formation histories. Those that are "dead" in that 
they have not experienced any recent star formation episodes 
comprise less than half of the total population. In many of 
the massive, bulge-dominated galaxies star formation occurs 
in the outer regions of the galaxy, where the light has sim- 
ply not been picked up by large scale redshift surveys that 
employ single fibre spectrographs. The fibres typically have 
diameters of 2-3 arcseconds and only sample light from the 
inner regions of the galaxies. 

The realization that many ellipticals do have extended 
star-forming disks first came from studie s of early-type 
galaxies at ultra-viole t wavelengths (e .g., lYi et alJ 120051 : 
iKauffmann et al.ll2007l : iFang et al.1 12012! ) . and these results 
are now bei ng confirmed by IF U studies of the nearby galaxy 
population IPerez et alJ (|2013). Our study reveals that the 
star formation in these disks did not occur continuously over 
a Hubble time, but was concentrated in a recent episode or 
burst. In future, large-scale Integ ral Field Unit (IFU) sur- 
veys of galaxies (see for example ISanchez" et alJ l2012h will 
shed more light on the nature of the star formation in the 
far outer reaches of nearby galaxies. 
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